Critical behavior of the frustrated antiferromagnetic six-state clock model on a 

triangular lattice 
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We study the anti-ferromagnetic six-state clock model with nearest neighbor interactions on a 
triangular lattice with extensive Monte-Carlo simulations. We find clear indications of two phase 
transitions at two different temperatures: Below Ti a chirality order sets in and by a thorough 
finite size scaling analysis of the specific heat and the chirality correlation length we show that this 
transition is in the Ising universality class (with a non- vanishing chirality order parameter below Tj). 
At Tkt{< Ti) the spin-spin correlation length as well as the spin susceptibility diverges according to 
a Kosterlitz-Thouless (KT) form and spin correlations decay algebraically below Tkt- We compare 
our results to recent x-ray diffraction experiments on the orientational ordering of CFaBr monolayers 
physisorbed on graphite. We argue that the six-state clock model describes the universal feature of 
the phase transition in the experimental system and that the orientational ordering belongs to the 
KT universality class. 



I. INTRODUCTION 

The study of frustrated two-dimensional spin models 
is motivated by their relevance for phase transitions in 
a wide range of physical systems. For instance, the 
fully frustrated XY (FFXY) model describes an ar- 
ray of Josephson junctions under an external magnetic 
field 1^,^ . Besides its experimental relevance the FFXY 
model has been the focus of a number of theoretical works 
since a delicate question arises here regarding the criti- 
cal behavior. The model has a continuous U{1) sym- 
metry and a discrete Ising, chiral or Z2 symmetry that 
can be broken at low temperatures through a Kosterlitz- 
Thouless (KT) type transition and an Ising-like transi- 
tion, respectively. 

Despite a lot of efforts [3-22], there is no consensus 
on the nature of the phase transitions and it is not 
clear whether the two transitions happen at two differ- 
ent temperatures or at a single one. Renormalization 
group (RG) studies on the FFXY model drew the con- 
clusion that the transitions occur at the same tempera- 
ture H]. Monte Carlo simulation studies on generalized 
FFXY models supported the single transition picture for 
the FFXY model Interestingly, it has been reported 
that the Z2 symmetry breaking transition may not be- 
long to the Ising universality class. Monte Carlo studies 
of the FFXY models on a triangular and a square lattice 
yielded the correlation length exponent ly — 0.83(4) (tri- 
angular lattice) and v — 0.85(3) (square lattice), which is 
inconsistent with the Ising value = 1 The specific 
heat appeared to follow a power-law scaling rather than 
the logarithmic scaling as one would expect for the Ising 
universality class. A non-Ising scaling behavior is also 
observed in the studies of the FFXY model via a Monte 
Carlo simulation and Monte Carlo transfer matrix cal- 



culations Jt],^- The same non-Ising critical behavior is 
also observed in the coupled XF- Ising model |5|,|9|,p^ . 

On the other hand, numerical evidence also in favor 
of two transitions at two different temperatures has been 
collected. Monte Carlo simulation studies of the frus- 
trated Coulomb gas system, which is supposed to be in 
the same universality class as the FFXY model, showed 
that the KT type transition temperature Tkt and the 
Ising-like transition temperature T/ are different with 
Tkt 



< Ti 
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Two transitions were also found in 
the FFXY model on a square lattice ||l^ and on a tri- 
angular lattice p^-p^ using Monte Carlo simulations. A 
careful analysis of the RG flow of the FFXY model also 
led to a conclusion of the double transition scenario . 
It is the general belief that the transitions at Tkt and 
at Ti belong to the KT universality class ||l^ and to the 
Ising universality class, respectively, if the transitions oc- 
cur at different temperatures. However, the critical expo- 
nents associated with the Z2 symmetry breaking, which 
are found by Monte Carlo simulations, turn out to be 
different from those of the Ising universality ||l^ , [l3| , p^ . 
They are rather close to those obtained in the XY-lsmg 
model Q. Although there is an argument that the ob- 
served non-Ising exponents are due to a screening effect 
hindering the asymptotic scaling behavior |l9|] , the con- 
troversy on the nature of the phase transition remains 
unsettled |^-||. 

Here we present a thorough numerical study of a re- 
lated model, the fully frustrated anti-ferromagnetic six- 
state clock model on a triangular lattice. We also re- 
port on an orientational ordering transition in a two- 
dimensional experimental system [ p3[ that, as we argue, 
is in the same universality class as the model we study 
numerically. The experimental system actually motivates 
(besides numerical simplicity) our restriction to the six 
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states of the spins rather than the continuum XY spins 
(and thus to a six-fold clock (Ce) symmetry rather than 
the U{1) symmetry). In the unfrustrated (ferromagnetic) 
case it has been shown, however, that the KT-behavior 
is stable with respect to a crystal field of six-fold sym- 
metry p4|-p6| . Hence we expect that our model displays 
the same critical behavior as the FFXY model. 

The paper is organized as follows: In Sec. |l^ we define 
the model and identify symmetry in the ground state. 
The corresponding order parameters are also defined. In 



Sec. in we explain briefly the Monte-Carlo procedure we 



used and determine the transition temperatures for the 
two phase transitions. Section IV is devoted to the classi- 



fication of the universal properties at the two transitions 
we find. In Sec. ^ we present experimental x-ray diffrac- 
tion results on the orientational ordering of CFaBr mono- 
layers physisorbed on graphite which we expect to be de- 
scribed by the KT-transition occurring in the model we 
studied numerically in the preceding sections. Section VI 
concludes the paper with a summary of our results. 



II. MODEL 

We investigate the phase transitions of the anti- 
ferromagnetic six-state clock model on a two-dimensional 
(2D) N — Lx X Ly triangular lattice (see Fig. |l]). The 
six-state clock spin S is a planar spin pointing toward 
discrete six directions; S — (cos 6*, sin ^^) with 



(n = 0,l,...,5). 



The interaction is given by the Hamiltonian 



n 
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FIG. 1. A triangular lattice of size 12 x 3. 
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The arrows 

represent the spin state in one of the ground states, Co- The 
lattice sites are labeled by the sublattice indices A, B, and C. 
Also shown are elementary up- and down- triangles denoted 
by Ai and Vii respectively. The ± signs denote the chirality 
of the triangles. 

increase by ^ along the up-triangles in the (anti-) clock- 
wise direction, and vice versa for down-triangles. 

There are two different symmetries between them; 
those in the same class of C or ^ are related by the rota- 
tion 9i ^ 6i + 2^ (m = 0, . . . , 5) and those in the others 
by the reflection 9i —9i. Therefore the model possesses 
the Ce (six-fold clock) symmetry and the Z2 (Ising) sym- 
metry. Note that the ¥YXY model on a triangular lattice 
has &U{1) and a Z2 symmetry Since the clock spin 
can have only six states, while the XY spin is a contin- 
uous one, the continuous U{1) symmetry of the ¥¥XY 
model is reduced to the discrete Ce symmetry. The na- 
ture of the Z2 symmetry in each model is identical. 

From the symmetry consideration, we expect that 
there are two types of phase transitions associated with 
the spontaneous breaking of the Cg symmetry and the 
Z2 symmetry. The chirality at each elementary triangle 
is defined as 



where the sum is over all nearest neighbor site pairs 
and J > is the antiferromagnetic coupling strength. 
The overall factor 2 is introduced for a computational 
convenience. 

The antiferromagnetic interaction on a triangular lat- 
tice induces a frustration. As a result the spins on each 
triangles should make an angle of ±120° with one another 
in the ground state. There are 12-fold degenerate ground 
states with three-sublattice structure. According to their 
chirality, we can categorize the ground states into C„ and 
■An (^1 = 0, 1, . . . , 5). In each ground state the spin con- 
figurations of the three sublattices A, B, C (Fig. |^) are 
given by 



27rn 
~6~ 



7ies — 



27r(ri± 2) 
"6 



^iec — 



27r(n±4) 
"6 



(3) 



for Cn (upper sign) and An (bottom sign). In other 
words, in the C{A) type ground state, the spin angles 



O + sin(0fe-0,) + sin(0,-0fe)] , 



where Aj and Vi are up- and down-triangles as depicted 
in Fig. |l|. Then the ground state C„ (An) has a checker- 
board pattern of the chirality with /ia = +1 (—1) and 
/ly = — 1 (+1). The staggered chirality 



1 

2N 



(4) 



plays a role of the order parameter for the Z2 symme- 
try breaking transition. The order parameter for the Cq 
symmetry breaking transition is the sublattice magneti- 
zation 



rriA = ^exp(i6',) 



(5) 



where the sum runs only over the sites in the A sublattice. 
niB and mc are defined analogously. 
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III. TRANSITION TEMPERATURES 

We performed Monte Carlo (MC) simulations on fi- 
nite N — L X {L/2) lattices with a sublattice updating 
scheme; one of the three sublattices is selected randomly 
and then all spins in the chosen sublattice are flipped 
according to the Metropolis rule. One Monte Carlo step 
corresponds to three sublattice updates. Various observ- 
ables are measured during the MC runs, such as the en- 
ergy, the chirality, and the sublattice magnetization, from 
which we can measure the averaged quantities and their 
fluctuations. In some cases, histograms are constructed 
from particularly long MC runs to obtain the observables 
as continuous functions of the temperature. 

The transition temperature is determined from the 
Binder parameter 



and Br, 



?/^2 \2 



(6) 



for the chirality and the magnetization, respectively. 
Here the angle bracket denotes a thermal average, which 
can be done by a time average over the MC runs. The 
Binder parameter is scale-independent at the critical 
points and approaches | (0) in the ordered (disordered) 
phase as L becomes larger. Hence the order-disorder 
transition temperature Tj related to the Z2 symmetry 
breaking is obtained from the crossing point in the plot 
of Bh versus T at different system sizes. From Fig. ^ (a), 
we estimate that 



Tj = 1.038 ±0.0005 . 

Bm also displays the crossing behavior (Fig. 
which we estimate that 

Tkt = 1.035 ±0.0005 



(7) 

I (b)), from 
(8) 



for the Cg symmetry breaking transition. 

The Binder parameters show different size dependence 
at low temperatures. For the staggered chirality Bh con- 
verges to 2/3 rapidly as L increases. However, B^ ap- 
pears to converge to values less than 2/3 at T < Tkt- 
It indicates that there is a quasi-long-range order in 
spins at T < Tkt- We confirm it from finite-size- 
scaling (FSS) behaviors of the magnetization. The sub- 
lattice magnetization shows the power law scaling be- 
havior, (|toa|) L~^, at T < Tkt with temperature 
dependent exponent x (see Fig. ||). 

The two transition temperatures lie very close to each 
other. Nevertheless the accuracy of the data is suffi- 
ciently high that we can exclude a possibility Tj = Tkt- 
Figure |^ demonstrates this via the correlation functions. 
We measure the correlations between the chirality of up- 
triangles and the magnetization of spins in the A sublat- 
tice displaced by a distance r in the vertical direction: 
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FIG. 2. Binder parameters for the staggered chirality in 

(a) and the magnetization in (b). The estimated transition 

temperatures are marked by broken lines. 
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FIG. 3. Sublattice magnetization (|?nA|) near the critical 
temperature Tkt- At temperatures T < Tkt it shows a 
power law behavior. The broken line has a slope —0.17. 
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Here, i+r denotes a site displaced vertically by a distance 
r from i, and denotes an up-triangle whose left corner 
is i (see Fig. |^) . In Fig. ^, we plot both correlation func- 
tions in the log-log scale at an intermediate temperature 
T = 1.037. Clearly one can see an upward curvature in 
the plot of Ch{r) for r <C Ly/2, which implies that the 
chirality order has already set in. On the other hand, 
there is a downward curvature in the plot of C„i(r) for 
r <C Ly/2 indicating that the magnetic order has not set 
in yet. Therefore we conclude that Tkt < 1.037 < Tj. 
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FIG. 4. Correlation functions of the staggered chirality, Ch-, 

and the magnetization, Cm, at temperature T = 1.037. 



IV. UNIVERSALITY CLASS 

With Tkt 7^ 7/, one expects that the both symme- 
try breakings take place independently. Then, the 
symmetry breaking transition should be in the Ising uni- 
versality class. The Ce symmetry is equivalent to that of 
the XY model perturbed by the (p = 6)-fold anisotropy 
field. The unperturbed XY model displays a KT tran- 
sition which separates a disordered high-temperature 
phase and a quasi-long range ordered low-temperature 
phase. A renormalization group (RG) study shows that 
an anisotropy field with p > 4 does not change the 
KT nature of the transition [Q. The same is true for 
the extreme case of the ferromagnetic six-state clock 
model ]2^ ]. So the Cg symmetry breaking phase tran- 
sition at r = Tkt is expected to belong to the KT uni- 
versality class. 

The interplay of the Ising type and the KT type or- 
dering has long been studied in the context of the fully 
frustrated XY models on a square or a triangular lattice. 
Frequently, it has been claimed [^|-^ that the Ising-like 
transition and the KT-type transition occur at the same 
temperature. On the other hand, more recent studies 
have reported the transitions to occur at two different 
temperatures JTT|- ^ , p^ , pOp^ ] . Surprisingly, some of re- 
cent high-precision Monte Carlo simulation studies sug- 
gested that the Z2 symmetry breaking transition does 
not belong to the Ising universality class in spite of the 
double transition 1 13 1^J2^. They estimated that v ~ 0.8 
for the correlation length exponent and ajv 0.46 with 
the specific exponent a, which are incompatible with the 
Ising values of = 1 and a = 0. The critical exponents 
are consistent with the values observed along the single 
transition line of the coupled XF-Ising model |]5|,p|,|l0| . 
Olsson [ p^ argued that the measured non-Ising critical 
exponents are artifacts of unusual finite size effects origi- 
nated from non-critical spin-wave fluctuations. However, 
apparently there is no general agreement yet on this sub- 



ject 

In our model, we observe the phase transitions occur- 
ring at two different temperatures. However, the contro- 
versy existing in the YYXY model tells us that it does 
not guarantee necessarily that the individual transitions 
will belong to the KT universality class and the Ising 
universality class, respectively. Therefore we perform a 
thorough finite-size-scaling (FSS) analysis to understand 
the nature of the phase transition. We believe that the 
antiferromagnetic six-state clock model can be used to 
resolve the existing controversy for the FFXF model. It 
has the proper symmetry properties as discussed previ- 
ously. And, from a practical point of view, larger system 
sizes are available since it is a discrete model, with which 
one can reduce the finite size effects. 



A. Z-i symmetry breaking 

We investigate the nature of the Z2 symmetry breaking 
transition by exploring the FSS property of the specific 
heat, defined by 



iV[(e 



(11) 



with e the energy per site. For a finite system (X — 
L X L/2) the specific heat has a peak c* ^ L"/" near 
the critical temperature, diverging with system size L — > 
00. Since the specific heat does not diverge at a KT 
type transition, the divergence is due to the Z2 symmetry 
breaking phase transition. The Ising universality class 
has a = 0, and the specific heat shows a logarithmic 
divergence, c* ^ InL. 

The specific heat is measured accurately using the 
standard Monte Carlo (MC) method combined with a 
single histogram method p7| . First we measure the spe- 
cific heat from the standard MC simulations at discrete 
temperature grid, which leads to a temperature Tq at 
which the specific heat is maximum approximately. Then 
the histogram 17o({S}), frequency of a spin state {S}, is 
constructed from MC runs at the temperature Tg over 
a time interval At. It enables us to obtain the average 
value of any observable Q as a continuous function of 
temperature T near Tq; 



E{s}^^o({S}) e- 



_(T~i-T-^)iS({S}) 



It is important to take At as large as possible to have 
good statistics. At 192 x 96 lattice, for example, we take 
At = 5 X 10^ which is ~ 2000 times of the relaxation 
time of energy-energy auto-correlation. The maximum 
system size we can simulate is 384 x 192. In that case we 
take At = 10^ which is ^ 1200 times of the relaxation 
time of energy-energy auto-correlation. Figure ^ shows 
the specific heat obtained from independent MC runs and 
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FIG. 5. Specific heat obtained from independent MC 
runs (symbol) and tlie liistogram metfiod (solid lines). To 
denotes the temperature at which the histogram is obtained. 
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FIG. 6. Peak values of the specific heat versus L in log-log 

scale (a) and semi- log scale (b). In (a) the solid line has a 

slope 0.5. 



from the histogram method. Their agreement near the 
peak is excellent, although they begin to deviate shghtly 
away from the peak. With this technique we can estimate 
c* very accurately. A statistical error for c* is estimated 
in the following way; we construct ten sub-histograms 
and as many values of c* from each of them. The error 
bar is taken as the half distance between the maximum 
and minimum values among them. 

The result is shown in Fig. ^. Interestingly, there is 
a crossover from a power-law scaling behavior c* ~ L^'^ 
(solid line in Fig. ^(a)) at L < to a logarithmic scaling 
behavior c* ~ InL (solid line in Fig. ^(b)) at L > 
with Lc — 96. It indicates that the specific heat has 
a logarithmic scaling behavior asymptotically, which is 
consistent with the Ising universality class. 

We also study the scaling of the chirality correla- 
tion length. We measured the correlation length using 



tion defined in Eq. yt 



■ Ch {t) with the correlation func- 
It is plotted in Fig. |^. We find 
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FIG. 7. Chirality correlation length. The values of the ex- 
ponent 1/ obtained from the fit at each lattice size are shown 
in the main panel and plotted in the inset. 



a strong dependence of ^/j on the system size. In addi- 
tion to the saturation of ^/i to 0{L) near and below the 
critical point, the dependence of the correlation length 
on L in the regime <C i is visible through a varying 
slope (cf.. Fig. 0). As a consequence, a fitting to a form 

— a{T — Ti)^^ at small system size will lead to a lower 
value of v than its asymptotic value. The fitted values 
of V at each lattice size are plotted in the inset of Fig. |^. 
It increases as the system size and approaches the Ising 
value ~ 1.0 in the infinite size limit. 

From the analysis of the specific heat and the chiral- 
ity correlation length, we conclude that the Z2 symmetry 
breaking transition indeed belongs to the Ising universal- 
ity class. Note that the effective specific heat exponent 
oijv 0.5 and the correlation length exponent v ~ 0.8 
appearing at small length scale arc comparable with the 
results of previous MC works claiming the non-Ising na- 
ture in the frustrated XY models ||,|,|l^,0|6|j2^j2| . 
Our results are fully consistent with Olsson's argument 
claiming that the non-Ising exponents reported by oth- 
ers are explained by a failure of finite-size-scaling at small 
length scale due to the screening length associated with 
the nearby KT transition jl^. The strong finite size ef- 
fect is overcome in this work since we could study large 
systems. We expect that one could observe the same 
crossover to the asymptotic Ising type scaling behavior 
in larger scale simulations in the frustrated XY systems. 



B. Ce symmetry breaking 

The KT transition is characterized by the essential sin- 
gularity of the correlation length and the susceptibility 
approaching the critical temperature. We measure the 
correlation length as 
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: 2.54e" 



and Xm = 4.20e" 



(13) 



from the magnetic correlation fmiction Cm in Eq. (10). 
One can also estimate the correlation length by fitting 
Cm{r) to a form with e"''/^™ . We tried with various func- 
tion forms but found that Eq. (|l^) yields the estimate 
most stable against statistical errors. The magnetic sus- 
ceptibility is measured from the fluctuation of the mag- 
netization, Xm = -^('71^). Figure ^ shows the plot of 
and Xm versus T — Tkt in log-log scale. Upward curva- 
ture in both plots indicates a divergence stronger than 
algebraic. We fit those quantities with the KT scaling 
form, U = ae''(^-^^^)"'' and x™ = a'e^'^^-^'^^r^ with 
y = fixed, whose results are drawn with solid lines 
in Fig. ||. It shows that the KT scaling behaviors are 
emerging only after L > 96. If y is also taken as a fitting 
parameter, one obtains y < 0.5 whose value depends on 
the fitting range; it increases as approaching the critical 
temperature. 

In the XY phase (T < Tkt) the spins have the quasi- 
long range order; the magnetization scales algebraically 
as (Iti^iI) with temperature dependent exponent 

X (see Fig. ^). The RG theory predicts that x = 1/8 at 
the KT transition point [Q. Our numerical data show 
that a;~0.17atr = Tkt- We think that the dis- 
crepancy is originated from a sensitive dependence of x 
on T. In summary, in spite of a quantitative disagree- 
ment, the Cg symmetry breaking phase transition is in 
qualitative agreement with the KT universality class; we 
observe the essential singularity in the correlation length 
and the susceptibility and quasi-long range order in the 
low temperature phase. 



V. EXPERIMENTAL REALIZATION: CF3BR ON 
GRAPHITE 

In this section we describe a possible experimental re- 
alization of the theoretical model we investigated above. 
In a recent work we reported results on x-ray powder 
diffraction study on a monolayer of halomethane CFaBr 
adsorbed on exfoliated graphite |Q. CFaBr is a pro- 
late molecule and has a dipole moment of about 0.5 D. 
The coverage p, temperature T phase diagram is rather 
complex [ p8|p9| ]. In we concentrated on a coverage 
which is representative of the extended monolayer regime 
in which the monolayer lattice is commensurate with the 
graphite lattice. This yields a 2 x 2 triangular lattice 
arrangement of the CFsBr molecules below a tempera- 
ture of 105 K which is the melting temperature of 
the commensurate layer. The inter-molecular distance is 
a — 4.92 A. Note that the lateral size of the graphite crys- 
tallites is only around 180 A, which confines any spatial 
correlation length to this value. 

An isolated CFaBr would prefer to lie fiat on the sub- 
strate, but the 2x2 mesh is too tight to accommodate 
the molecules in this orientation. Therefore the individ- 
ual CFsBr molecules stand on the substrate, presumably 
with the F3 tripod down, with a maximum tilt angles of 
the molecular axis up to 30° with respect to the substrate 
normal due to steric repulsion. A tilt leads to a non-zero 
in-plane component of the dipole moment. We regard 
this component as planar pseudospin Si = (cos 6*^, sin 0^) 
with the azimuthal angle 9i of the molecule. In this sense 
the 2x2 state is disordered with a zero time average of 
every Si, and is stabilized at higher temperatures by a 
libration and/or a precession of the molecular axis about 
the substrate normal. 

As the temperature is decreased additional features de- 
velop in the diffraction pattern which finally, below 40 K, 
can be identified p3| as Bragg peaks (with a finite width 
of around (180 A)"^^ due to the lateral size of the crys- 
tallites) indicating an orientational order in the dipole 
moments identical to the one depicted in Fig. |l|. The 
temperature dependence of the correlation length ^, de- 
termined from the intrinsic width of this peak, can be 
fitted with the KT-expression 



^ = Acxp{BiT/TKT-iy'^^) 



(14) 



to the data for T > 40 K (see Fig. 3 of Ref. [g3|). The 
fit parameters are A = 9 ± 2 A, B = 1.5 ± 0.4, Tkt = 
30 ± 3 K. Note that the value of A is reasonably close to 
the lattice parameter of the 2D mesh. Thus ^ is expected 
to diverge at a KT-critical temperature Tkt of about 30 
K, but the growth of the correlated regions is interrupted 
when ^ reaches the size of the graphite crystallites. This 
happens at about 40 K. 

We think that the model (||) is a good description of the 
orientational ordering process described in this physical 
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FIG. 9. Comparison of the magnetic correlation length 
(scaled by the lattice constant) from the experimental data 
(from Fig. 3 of Ref. |^) and the Monte Carlo simulations 
(from Fig. ^). The reduced temperature {T/Tkt — 1) for 
the simulation data are rescaled by a factor of 18 in order to 
achieve an acceptable data collapse for linear lattice sizes L 
between 96 and 192. 



system: Clearly the pseudospin correlations are bound to 
a plane, thus the system is 2D with respect to the relevant 
degrees of freedom at the phase transition. Moreover, as 
mentioned before, below 105 K the CFsBr molecules are 
arranged in a triangular lattice. The pseudospin repre- 
senting the CFsBr dipole moment is presumably not a 
strictly isotropic planar rotator but experiences a crys- 
tal field from the graphite substrate which breaks the 
continuous azimuthal symmetry into the six fold symme- 
try of the monolayer. This is reflected by the six-state 
clock variables of model (||). The ordered structure of 
the CFaBr dipoles is antiferroelectric, which is taken into 
account by the antiferromagnetic couplings between the 
pseudospins in (H). Finally, the character of the rele- 
vant orientation dependent interactions of CFaBr is not 
known, but a comparison of monolayers of several po- 
lar methane derivatives ||2^ shows that fully halogenated 
molecules including CFaBr with small dipole moments 
around 0.5 D have structures different from partially 
halogenated molecules such as CH3CI with strong dipole 
moments around 1.7 D. This suggests that for CFsBr the 
short range anisotropic part of the intcrmolccular van der 
Waals force and hard-core repulsion are more important 
than the medium range dipole-dipole interaction. Thus 
the interactions can be assumed to be short ranged. Thus 
one expects that model (^) and the physical system dis- 
cussed here are in the same universality class. 

In Fig. ^ a comparison of the magnetic correlation 
length from the experiment p3| and the simulation (from 
Fig. ^) is shown. Since the factor B in the KT-form ([T^ ) 
of the correlation length is a non-universal number the re- 
duced temperature {T/Tkt — 1) has to be rescaled by an 
appropriate factor in order to achieve an acceptable data 



collapse. The rescaling factor turns out to be quite large, 
namely 18, which is not unusual for microscopically dif- 
ferent systems in the KT universality class (see e.g., PH). 
Note that the finite linear size of the crystallites plays a 
similar role as the finite lattice sizes in the simulations 
and sets the saturation value for the correlation length 
(divided by the lattice constant for the triangular lattice 
of the CFsBr molecules, which is a = 4.92 A). The nice 
collapse of the experiment and simulation data supports 
our claim that the two physical systems are in the same 
universality class and that the orientational ordering of 
CFsBr molecules on graphite is in the KT universality 
class. 



VI. SUMMARY 

To summarize we have studied the phase transitions 
in the anti-ferromagnetic six-state clock model on a tri- 
angular lattice, which is fully frustrated. As a result 
the ground states have a Cg (six-state clock) symmetry 
and a Z2 (Ising) symmetry. Through extensive Monte 
Carlo simulations we found that the model undergoes 
a Kosterlitz-Thouless transition at Tkt and an Ising 
transition at T/. The two transitions correspond to the 
Cq and the Z2 symmetry breaking transition, respec- 
tively. High-precision Monte Carlo data indicate that 
the two transitions take place at different temperatures, 
Tkt < Ti (Eqs. (0) and (|)). This has been checked ex- 
plicitly by analyzing the behavior of the spin and the chi- 
rality correlation function at temperatures between Tkt 
and Ti (Fig. ^). Furthermore, we have shown that the Z2 
symmetry breaking transition belongs to the Ising univer- 
sality class. For small system sizes, the scaling property 
of the specific heat and the correlation length deviates 
apparently from the Ising universality class. However, 
simulation results for larger system sizes indicate that 
the model belongs asymptotically to the Ising universal- 
ity class (Figs. ||and0). As for the transition at Tkt, we 
have found that the magnetization correlation length and 
the susceptibility diverges at T = Tkt according to the 
KT scaling form (Fig. ^). We have also found that the 
spins have a quasi long range order below Tkt (Fig. |^). 
Combining these, we conclude that the transition at Tkt 
belongs to the KT universality class. 

Our model is a variant of the fully frustrated XY mod- 
els where the KT type ordering and the Ising type order- 
ing interplay interestingly. Our numerical results support 
a scenario that there are two separate phase transitions 
with Tkt 7^ Tj; one at Tkt in the KT universality class 
and the other at T/ in the Ising universality class. Our 
results are consistent with those in Ref. very well. 

We have proposed that our theoretical model describes 
the orientational ordering transition of CFsBr molecules 
on graphite since the model has the same symmetry as 
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the experimental system. We argue that the orienta- 
tional ordering transition belongs to the KT universality 
class by comparing the magnetic correlation length 
obtained from the experiment p3t with the correlation 
length obtained numerically in the six-state clock model. 
With a suitable rescaling of parameters, we have shown 
that the correlation lengths in both systems have the 
same scaling behavior (Figj9|). It gives more evidence 
that the orientational ordering transition is indeed the 
KT transition. 
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